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First Signs for String Breaking in Two-Flavor QCD 

Carleton DeTar ^, Urs Heller ^, Pierre Lacock*^* 

^Physics Department, University of Utah, Salt Lake City, UT 84112-0830, USA 
''SCRI, Florida State University, Tallahassee, FL 32306-4130, USA 

We have been examining the phenomenon of string breaking in QCD with two flavors of dynamical staggered 
quarks. We construct a transfer matrix from a combination of "string" and "two-meson" channels. Preliminary 
results with low statistics show the expected signs of string breaking. 



1. INTRODUCTION 

During the past few years there has been re- 
newed interest in the important phenomenon of 
string breaking (SB), which is predicted by QCD, 
but which lattice QCD simulations for a long time 
have failed to show conclusively. 

String breaking is observed as a leveling off of 
the static quark-antiquark potential at large sepa- 
ration. The potential is the separation-dependent 
ground-state eigenvalue of the QCD hamiltonian 
in the presence of the static quark-antiquark pair. 
Traditionally, this eigenvalue was sought in the 
Wilson loop observable, which is the expectation 
value of the transfer matrix on a "string" state 
with the string of color flux connecting the static 
quark and antiquark. Experience has shown, not 
surprisingly, that the string state (S) is a very 
poor variational ansatz for a state that looks more 
like two static-light mesons. Including an admix- 
ture of a two-meson component (M) should help 



In the following we report on preliminary re- 
sults that have been obtained using the rele- 
vant string-string and string-meson operators for 
QCD. 

2. STUDYING STRING BREAKING ON 
THE LATTICE 

Our conventional Wilson loop is computed with 
APE smearing of the space-like gauge links. In 
hamiltonian language the expectation value of 



this operator is the transfer matrix G'5s(-R, T) 
between an initial and final state S consist- 
ing of a static quark-antiquark pair separated 
by a fat string of color flux. We enlarge the 
space by including a meson-antimeson state M 
with an extra light quark in the vicinity of the 
static antiquark and an extra light antiquark in 
the vicinity of the static quark. Thus we also 
compute the additional transfer matrix elements 
Gmm{R.T), Gms{R,T) and Gsm{R,T). They 
are diagrammed in Figures 1 and 2. 



Figure 1. The static-light meson-antimeson pair 
contribution to the full QCD propagator. The 
wiggly lines denote the light quark propagator. 
Shown are the 'direct' and 'exchange' terms re- 
spectively. 

In principle, both channels couple to a common 
set of eigenvalues. At large T we expect to reach 
the ground state, defined by the largest general- 
ized eigenvalue 

G(i?, T + l)u{R, T) ^ \{R, T)G{R, T)u{R, T) (1) 
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The potential is then given by V{R) = 



2 



-liniT^oolog|A(i?,r)|. The vector u{R,T) de- 
fines the variationally optimum admixture of S 
and M with the largest overlap with the ground 
state. 

A key unitarity condition is that the eigen- 
values and eigenstates approach a constant with 
increasing T. This condition together with a 
demonstration of a smooth transition with in- 
creasing R between a string-dominated state and 
two-meson-dominated state constitutes a true 
test of string breaking and should distinguish a 
quenched calculation from a proper calculation 
with dynamical quarks. 

3. NUMERICAL METHOD 

To maximize statistics we generate "all-to-all" 
propagators for the light quark, using a Gaussian 
random source method. Results reported here are 
based on 15 such sources per gauge configuration, 
but we plan to increase this number |9|. 



Figure 2. The string- meson transition correlator 
GsM (and its hermitian conjugate Gms)- The 
wiggly line again denotes the light quark propa- 
gator. 

When constructing operators involving both 
staggered fermions and gauge links, one has to 
pay careful attention to staggered fermion phases. 
A consistent treatment results from interpreting 
the static quark as an infinitely massive staggered 
quark. For example, a hopping parameter expan- 
sion around an on-axis R x T rectangular path 
gives, in addition to the Wilson-loop gauge-link 
product, a net phase factor (—1)-'^-^ x (—1)-'^+^, 
independent of the Dirac phase conventions. We 
use a similar construction to get the phases for the 



nonclosed gauge link products in the diagrams of 
Figures 1 and 2. In that case the Dirac phase 
convention for the gauge link products must be 
consistent with that of the light quark. A pecu- 
liar consequence of this construction is that the 
transition matrix elements must vanish for off- 
axis displacements R that have more than one 
odd Cartesian-displacement component. 

4. LATTICE SIMULATION 

We use a set of stored configurations of size 
20^ X 24 at /3=5.415 generated with two flavors 
of dynamical staggered fermions of mass ma = 
0.0125. This set gives a ratio mT^/mp k, 0.358 and 
has a lattice spacing of about 0.17 fm, which gives 
a spatial lattice size of w 3.4 fm, and a temporal 
size of « 4 fm. 

For the light fermions in the static-light mesons 
we use the same parameters as for the dynamical 
fermions. With the choice of 15 random sources 
we have currently analyzed about 70 out of the 
200 archived configurations. 
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Figure 3. The behavior of the two- meson to two- 
meson matrix element Gmm as a function of the 
spatial separation R. Results were obtained from 
fits to the data as a function of T for fixed R. 

5. RESULTS 

In Figure 3 we show the results for the two- 
meson to two- meson transition Gmm- For short 
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distances the operator behaves hke the Wilson 
loop, with a linearly rising value, while for dis- 
tances greater than 3 in lattice units it starts to 
level off. The dashed line is twice the mass of the 
heavy-light meson. 

The results for the transition matrix clement 
Gms shown in Figure 4 are similar within large 
errors, with the energy leveling off at Ra w 5 in 
lattice units. 

In Figure 5 we show the combined results — 
including those for the pure Wilson loop correla- 
tor. The string-like behavior of the Gmm and, 
especially, of the Gms correlators at short dis- 
tances is clearly visible. 
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Figure 4. The same as in Fig. 2, but for the 
Gms operator. 



6. CONCLUSIONS 

Our present results with low statistics show 
a behavior expected with string breaking. The 
cross-over region is found to be at a distance of 
Ra = 5 — 6, or about 0.8 — 1.1 fm. However, 
to demonstrate string breaking convincingly, one 
must find a smooth transition with increasing R 
between a string-dominated state and two-meson- 
dominated state and demonstrate that the eigen- 
values and eigenvectors of the transfer matrix ap- 
proach a constant at large T. Work in this direc- 



tion is currently underway. 
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Figure 5. Combining the results for the Gmm 
(□), Gms (o) and Wilson loop {Gss) (x) corre- 
lators. 
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